Random bond XXZ chains with modulated couplings 
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The magnetization behavior of g-periodic antiferromagnetic spin 1 /2 Heisenberg chains under uni- 
form magnetic fields is investigated in a background of disorder exchange distributions. By means 
of both real space decimation procedures and numerical diagonalizations in XX chains, it is found 
that for binary disorder the magnetization exhibits wide plateaux at values of 1 + 2(p — where 
p is the disorder strength. In contrast, no spin gaps are observed in the presence of continuous 
exchange distributions. We also study the magnetic susceptibility at low magnetic fields. 
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PACS numbers: 75.10.Jm, 75.10.Nr, 75.60.Ej. 

The study of low dimensional antiferromagnets has 
received much renewed attention largely owing to the 
synthesis of ladder materials [|J. One particular issue 
that captured both experimental and theoretical efforts 
in the past few years is the appearance of magnetization 
plateaux, i.e. massive spin excitations in the magnetiza- 
tion curve. In general, the latter are quite robust and for 
pure spin systems appear at rational magnetization val- 
ues 1^-^]. More recently, some experiments have indeed 
confirmed these theoretical predictions in a few particu- 
lar cases Q but some issues yet remain unresolved . 

In order to make closer contact with experiments, one 
has to take into account the disorder that is almost in- 
evitably present due to lattice imperfections and mag- 
netic doping. A relevant question related to the appear- 
ance of magnetization plateaux is whether they are ro- 
bust in the presence of quenched disorder. As a first step 
in this direction, in this work we analyze the effect of a 
disordered distribution of exchange couplings with a peri- 
odically modulated mean on the magnetic behavior of an 
XXZ antiferromagnetic chain. Both even j?| and odd 
H modulations are known to exist and are ultimately 
responsible for the structure of the magnetization curve 
[ 4pp[| . y-merized XX chains have also been studied in 
[ HiT using a Jordan- Wigner transformation. 

By means of a decimation procedure similar to that 
used in jll],[l2| we argue that plateaux in the magneti- 
zation curve appear at specific magnetization values m 
which depend both on the couplings periodicity q and 
the strength of the disorder p. Hence disorder, instead of 
removing completely the plateau structure, shifts the po- 
sition of certain plateaux in a precise way which depends 
on the disorder strength. Surprisingly, as we shall see 
from our numerical evidence, the plateaux predicted via 
this simple argument are indeed present. Moreover, they 
are rather wide and therefore could be eventually observ- 
able in low temperature experiments which in turn would 
allow for a precise determination of the disorder degree. 

By extending the methods of we also investigate 



the characteristics of the magnetic susceptibility at low 
fields. Its behavior shows an interesting even-odd effect, 
similar to that found in the study of disordered XX TV- 
leg ladders |l4[ . In fact, for q odd we find the same kind 
of diver gen ce as for the homogeneously disordered case 
(? = 1)@, namely 
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whereas the even q modulations yield a generic non- 
universal power law behavior as in [p~5|Jlq] for q = 2, 
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where we give an analytic expression for a for arbitrary 
q. In principle, these results should emerge from experi- 
mental susceptibility measurements in disordered low di- 
mensional compounds. 

In what follows we focus attention on the occurrence 
conditions of zero temperature plateaux in a random q- 
merized XXZ spin-1/2 chain whose Hamiltonian is 

H = E J * ( S ? S Ui + S ? S ?+i + AS * S Ui ) > (3) 



where Ji are randomly distributed bonds. Specifically, 
let us consider a binary distribution of strength p (p — 
corresponds to the pure g-merized case while p = 1 
corresponds to the uniform chain), 



P{Ji) = P 5{Ji - J 1 ) + (1 -p) 6(Ji - Jo - ji J) 



(4) 



where 7$ = 7, (—7) if i = qn, (i ^ qn), along with 



a Gaussian disorder P(Ji) oc exp - 



and log- 



normal distribution given in terms of Wi — ln(Ji) and 



P(Wi) oc exp- 
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All these distributions, taken 



with same mean and variance, are built to enforce q- 
merization, whose value is measured on average by 7J. 
In what follows we assume that J' is the smallest coupling 
and consider < 7 < Jq/J. 
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a. Decimation procedure — Here we follow the argu- 
ments used by Fisher in [fL2| . The procedure is roughly 
to decimate the strongest bonds up to an energy scale 
given by the temperature. The remaining spins can be 
considered as free and each of them will then give rise to 
a Curie behavior in the magnetic susceptibility. 

In our problem (which is at T = 0) the energy scale is 
provided by the magnetic field, and in order to compute 
the magnetization, decimation has to be stopped at an 
energy scale of the order of the magnetic field. We as- 
sume that all spins coupled by bonds stronger than the 
magnetic field form singlets and do not contribute to the 
magnetization, whereas spins coupled by weaker bonds 
are completely polarized. The magnetization is thus pro- 
portional to the fraction of remaining spins at the step 
where we stop decimation. Our argument happens to ap- 
ply well to the binary distribution, provided the energy 
scales of the involved exchanges are well separated. 

b. Plateaux in q-merized chains — Let us first con- 
sider the case q = 2 and assume we start at high enough 
magnetic field, such that all spins are polarized (satura- 
tion, m = 1) and begin decreasing the magnetic field. 
The magnetization stays constant for a while, then de- 
creases abruptly at h ~ Jo + 7 J and after that a plateau 
occurs at m = p. This can be easily understood: at 
h ~ Jo + 7 J we can decimate all the strongest bonds 
Jo + 7 J (the corresponding spin pairs form singlets and 
do not contribute) , and the number of remaining (com- 
pletely polarized) spins is N — 2 x (1 — p)N/2 — pN . 
Here, the factor of 2 comes from the removal of two spins 
each time we decimate a bond. Hence, the first plateau 
occurs at to — p. The appearance of this spin gap is due 
to the fact that the remaining strongest bonds have val- 
ues Jo — 7 J p7|] , and all spins left from the first step of 
decimation remain polarized (and the magnetization con- 
stant), until the magnetic field decreases to h ~ J — 7 J. 
At this point the magnetization again decreases abruptly 
and a second plateau occurs. The abrupt change corre- 
sponds to the decimation of bonds Jo — 7 J which leaves 
us with N - 2 x (1 -p)^ - 2 X (1 - p)§ p 2 completely 
polarized bonds. The plateau occurs then at magneti- 
zation to = p — p 2 + p 3 . The term (1 — p)^p 2 comes 
from the bonds Jo — 7 J which, having a J' bond at each 
side, were not decimated in the first step and thus is the 
number of bonds actually decimated at the second step. 
Evidently, for q > 2 we can follow the same reasoning. 
Thus, the number of spins which yield finite contribu- 
tions to the total magnetization at h ~ Jo + 7 J is simply 
N — 2 x (1 — p)N/q. Hence, we find the first plateau at 

m = l + -(p-l). (5) 

q 

Notice that this result locates correctly the spin gaps ap- 
pearing in a pure g-merized chain (p = 0). In this sense, 
Eq. (||) provides an extension of this latter case || in the 
presence of binary disorder. Since the decimation proce- 
dure applies for generic XX Z chains p^ [, we conclude 
that the emergence of the plateaux predicted in d|) is a 



generic feature, at least with the anisotropy parameter 
A I < 1. It is straightforward to generalize this analysis 
to the case of an arbitrary but discrete probability dis- 
tribution. Given a finite difference between the highest 
values of the couplings in the inequivalent sites, one can 
predict the presence and position of the plateaux. 

To enable an independent check of these assertions, we 
turn to a numerical diagonalization of the Hamiltonian 
contenting ourselves with the analysis of the subcase 
A = 0. This allows us to explore rather long chains using 
a fair number of disorder realizations. In Figs. 1(a), 1(b) 
and 1(c) we show respectively the whole magnetization 
curves obtained for q — 2, 3 and 4 after averaging over 
100 samples of L — 5 x 10 4 sites under the exchange dis- 
order (H). It can be readily verified that this set of robust 
plateaux emerges quite precisely at the critical magneti- 
zations given by Eq. (j^) . The secondary plateaux, though 
narrower, are still visible in Fig. (1). 

It is important to stress that the derivation of our re- 
sults for the quantization conditions derived above rely 
strongly on the discreteness of the probability distribu- 
tion and would not be applicable to an arbitrary contin- 
uous exchange disorder. In fact, for the Gaussian case 
referred to above it turns out that no traces of plateaux 
can be observed. Furthermore, the magnetic suscepti- 
bility in the Gaussian case only vanishes asymptotically 
when approaching the saturation regime, as can be seen 
in Fig. 2 for a variety of coupling periodicities. Here, the 
sampling was improved up to 5 x 10 4 realizations though 
the length of the chain was reduced to L = 1500 sites, as 
the CPU time per spectrum grows as L 2 . 

^From the numerical curves, it appears that the usual 
DN-PT |l^] transition is smoothed in the presence of 
disorder, both in the cases of binary and Gaussian dis- 
tribution. Using the exact results of [|l9| for a fam- 
ily of Poissonian distributions, one sees that the behav- 
ior of the magnetization close to saturation has a non- 
universal exponential decrease. On the other hand, this 
non-universality is reflected for the binary case, in the 
fact that saturation occurs with an upper bound given 
by 2J max . Since the universal DN-PT transition is de- 
stroyed near saturation, we expect that the same will 
occur in the vicinity of a non-trivial plateau. In fact, 
this is noticeable in the numerical data. 

c. Susceptibility at low magnetic fields — For homoge- 
neously disordered chains (i.e. q = 1), one can use the 
decimation procedure of |1^] along with the universal- 
ity of the fixed point, to show that either for discrete or 
continuous distributions the low field magnetic suscepti- 
bility behaves according to Eq. ([!]). Following a simple 
argument based on random walk motion used in it 
can be readily shown that for A = (or XX chains), 
these arguments can be extended to the case of q odd 
giving the same singularities. In fact, these expectations 
can be compared to the numerical results obtained for q 
odd with both Gaussian and binary disorders, as shown 
in Figs. 2 and 3 respectively. In particular, we direct 
the reader to the semi-log insets of Figs. 2(b), 2(d), 3(b) 
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and 3(d) which evidently follow the universal singularity 
referred to in Eq. (Q). The numerical results for the log- 
normal distribution lead to the same qualitative behavior 
obtained for the Gaussian case. 

For q even, for which there is a plateau at m = in 
the pure case, the situation is more subtle. Using the 
notation of Jl3|] , for XX chains we can again define a 
random walk of the variable Ui — ln(Aj) between the 
boundaries ln(V 2 /E) and hi(E), now with a driving force 
F and diffusion coefficient D given by 

2 

F = - < In Ji =qn - In Jf^ qn > , (6) 

D = - [var 2 (ln jf =qn ) + (q - l)var 2 (ln J^j] . (7) 

By means of the method given in |i~3[ ] for the undriven 
random walk, one can approximate the problem to a dis- 
crete time diffusion problem with an absorbing and re- 
flecting wall. One can then show that the average num- 
ber of bonds for a cycle to be completed goes now like 
n ~ e a A "/ 2 , which gives the asymptotic behavior for the 
magnetic susceptibility as in (0). The non- universal ex- 
ponent a turns out to depend on the distribution param- 
eters (^) and (^) namely, a = 2F/D. Also, this exponent 
coincides with the results obtained in |l6[] using decima- 
tion and other methods for the q = 2, XXZ chain. 

Once more, our numerical results for the binary and 
Gaussian coupling distributions considered above lend 
further support to this non-universal picture of even 
modulations. Specifically, there are in fact situations 
for which the low field susceptibility can either diverge 
(a < 1) as displayed in Figs. 2(a), 2(c), 3(c) or collapse 
(a > 1) as shown in Figs. 3(a). Moreover, we checked 
that our a exponents can fit reasonably the numerical be- 
havior obtained in this regime, as indicated by the log-log 
insets of Figs. 2(a) and 2(c). 

To summarize, we have studied the effect of disor- 
der on the plateaux structure in g-merized XXZ chains. 
By means of a simple real space decimation procedure 
we could account for a non-trivial phenomenon namely, 
the shift in the magnetization values for which certain 
plateaux emerge, as compared to the pure system. This 
was tested by numerical diagonalizations of large XX 
chains finding a remarkable agreement with Eq. (H) . We 
have also analyzed the behavior of the low magnetic field 
susceptibility which exhibits a clear g-odd (even) loga- 
rithmic (power law) behavior. Our theoretical predic- 
tions could be experimentally checked on dimerized com- 
pounds such as CuGeOs doped with Si [^0| under mag- 
netic fields. Also trimerized compounds exist in nature 
H and it would be interesting to see if they can be doped. 
We trust this work will convey an interesting motivation 
for further experimental studies. 
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FIG. 1. Magnetization curves of modulated XX spin 
chains with q = 2 (a), 3 (b) and 4 (c), immersed in 
disordered binary backgrounds of strength p. Solid lines 
represent averages over 100 samples with 5 x 10 4 sites, 
J' /Jo = 0.2, 7J = 0.5 and p = 0.2, 0.4, 0.6, 0.8 in ascending 
order. The left and rightmost dotted lines denote respectively 
the pure uniform and pure modulated cases p = 1 and p = . 
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FIG. 2. Magnetic susceptibility of g-merized XX chains 
after averaging over 5 x 10 4 samples with 1500 sites, using 
q-periodic Gaussian exchange distributions for q — 2 (a), 3 
(b), 4 (c), 5 (d) and strength p — 0.4. The insets show respec- 
tively the susceptibility behavior at low magnetic fields which 
follows closely the regimes predicted by Eqs. (1) (q odd), and 
(2) (q even), in the text. 
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FIG. 3. Same as Fig. 2, but averaging over a binary ex- 
change disorder of strength p — 0.4. The small field sus- 
ceptibility behavior displayed in the insets reflect the typi- 
cal singularity of odd periodicities [ q = 3 (b) and 5 (d) ] , 
whereas even modulated distributions [ q = 2 (a) and 4 (c) ] , 
are non-univesal in this regime. 
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